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Abstract 


A  simple  lightly  damped  second-order  system  is 
augmented  with  an  optimal  controller  using  linear  quadra¬ 
tic  regulator  theory.  However,  instead  of  using  con¬ 
ventional  optimal  control  design  methods,  the  controller 
is  designed  with  frequency  weighted  cost  functionals. 

This  report  investigates  the  effect  of  frequency  shaping 
on  the  frequency  and  time  response  characteris tics  of 
the  closed-loop  system. 

The  effects  of  frequency  shaping  on  the  state  and 
control  penalty  matrices  are  analyzed  separately.  Each 
test  case  yields  improved  damping  ratios  compared  with 
the  baseline  open-loop  system.  Each  system,  including 
the  baseline  model,  has  a  high  frequency  asymptote  slope 
of  -40  dB/decade.  There  is  very  little  difference  in  the 
high-frequency  response  characteristics  of  the  frequency 
shaped  augmented  systems  and  the  open-loop  system. 

The  most  important  result  is  that  the  low  frequency 
magnitude  response  can  be  reduced  by  using  frequency 
shapings.  This  is  particularly  useful  in  attenuating 
low  frequency  system  modes  as  well  as  low  frequency  noise. 
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OPTIMAL  CONTROL  USING 
FREQUENCY  WEIGHTED  COST  FUNCTIONALS 

I  Introduction 

Control  system  design  has  typically  involved  either 
frequency  domain  compensation  or  state  space  time  domain 
techniques.  Using  frequency  domain  analysis  techniques, 
control  systems  designers  have  developed  an  understanding 
and  an  intuition  for  what  type  of  frequency  response 
characteristics  a  system  should  have.  The  design  tech¬ 
nique  is  then  to  add  compensation  in  such  a  manner  so  as 
to  obtain  this  desired  frequency  response. 

An  advantage  of  frequency  domain  analysis  is  that 
there  exists  a  great  body  of  knowledge  concerning  this 
method.  For  example,  the  use  of  washout  circuits  in 
aircraft  yaw  dampers  is  very  common,  very  effective,  and 
done  purely  in  the  frequency  domain. 

A  significant  disadvantage  to  frequency  domain 
analysis,  which  is  discussed  in  Ch.  2,  is  that  it  is  dif¬ 
ficult  to  apply  the  method  to  multiple-input-multiple- 
output  (MIMO)  systems. 

State  space  techniques  differ  in  that  they  deal 
directly  with  the  time  domain  characteristics  of  the 
system.  A  common  state  space  technique  is  the  use  of 
linear  quadratic  optimal  control.  The  designer  defines 
a  quadratic  cost  functional  and  then  by  minimizing  the 
cost  functional  determines  the  optimal  control.  One  main 
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advantage  to  the  state  space  method,  discussed  in  Ch.  2, 
is  the  ease  at  which  it  can  be  applied  to  MIMO  systems 
through  the  use  of  vector  and  matrix  representations. 
Optimal  control  design  generally  requires  choosing  con¬ 
stant  weighting  matrices  for  the  cost  functional. 

There  has  been  much  research  on  determining  the 
control  of  a  system  using  both  the  frequency  domain  and 
time  domain  systems.  However,  it  is  not  clear  what 
advantages  might  be  gained  by  using  some  combination 
of  the  two  methods. 

Background 

Recent  work  by  N.  K.  Gupta  (Refs  6,7)  has  shown 
that  optimal  control  theory  can  also  be  applied  to  cost 
functionals  involving  frequency  dependent  weighting 
matrices.  Gupta  suggests  that  this  method  is  very  helpful 
in  controlling  spillover,  reducing  high  frequency  response, 
and  improving  disturbance  rejection  (Ref  6) .  It  is  cur¬ 
rently  known  how  to  formulate  this  problem,  but  Gupta 
has  not  shown  or  suggested  any  relationship  between  the 
cost  functionals  and  the  frequency  response  of  the  closed- 
loop  system. 

Problem 

The  problem  is  to  investigate  the  relationship 
between  these  frequency  weighted  cost  functionals  and  the 
conventional  frequency  response  analysis  techniques. 
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Specifically,  frequency  response  characteristics  such  as 
magnitude  overshoot,  high  frequency  asymptote  slope,  break 
frequency,  and  damping  are  examined.  Additionally,  time 
response  characteristics  of  the  systems  are  examined. 

Assumptions 

It  is  assumed  that  all  system  models  are  linear 
deterministic,  that  is,  that  the  systems  are  noise-free. 

General  Approach 

Using  a  simple  plant,  several  frequency  shaped  cost 
functionals  proposed  by  Gupta  are  analyzed. 

The  closed-loop  response  of  a  second  order  baseline 
system  with  no  controller  is  analyzed.  Then  the  baseline 
system  is  augmented  with  a  baseline  controller.  Subsequent 
cases  examine  the  closed-loop  response  of  the  baseline 
system  augmented  with  Gupta's  suggested  frequency  weighted 
controllers.  The  frequency  weighted  controllers  are  com¬ 
pared  with  the  baseline  controllers  in  terms  of  frequency 
and  time  response.  The  effect  of  frequency  shapings  on 
both  the  state  and  the  control  penalty  matrices  is 
examined . 
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II  Theoretical  Development 


Many  feedback  control  systems  can  be  designed  using 
the  linear  time  invariant  system  model 

x(t)  =  Ax  (t)  +  Bja(t)  x(o)  =  x^  (o<t<t£)  (2-1) 

y(t)  =  Cx ( t)  +  Du ( t)  (2-2) 

where  x  is  an  n-dimens ional  state  vector,  u  is  an  Tri¬ 
dimensional  control  vector,  y  is  a  p-dimensional  output 
vector,  A  is  an  n  x  n  plant  matrix,  B  is  an  n  x  m  con¬ 
trol  matrix,  c  is  a  p  x  n  output  matrix,  and  D  is  a 
p  x  m  feeu -forward  matrix. 

The  control,  u,  can  be  defined  to  be  the  product  of 
a  constant  feedback  matrix,  K,  and  the  state  vector,  x(t) , 
in  the  following  manner: 

u(t)  =  -Kx ( t )  (2-3) 

There  are  various  methods  of  determining  the  feedback  gain 
matrix,  K.  Pole  placement  techniques,  entire  eigen- 
structure  assignment,  steady  state  solution  of  appropriate 
Riccati  equations  via  linear  quadratic  methodology,  and 
even  arbitrary  selection  are  some  of  these  methods.  The 
value  of  K  is  manipulated  until  the  desired  closed-loop 
system  dynamics  are  achieved.  Fig.  2-1  shows  a  block 
diagram  of  the  closed-loop  system  where  r  is  an  arbitrary 
input.  It  r  is  assumed  to  be  a  constant  value  other  than 
zero  and  the  feedback  u(t)  =  -K  x(t)  +  r(t)  is  used, 
the  state  equation  for  the  closed-loop  system  becomes 
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x(t)  =  Ax(t)  +  Bu(t)  =  (A  -  BK)x(t)  +  Br  (t)  (2-4) 

where  u(t)  =  r (t)  -  Kx(t) .  The  poles  of  the  system  trans¬ 
fer  function  between  x(t)  and  u(t)  are  the  eigenvalues  of 
(A  -  B  K)  . 

Optimal  Control  Methods 

Using  linear  optimal  theory,  the  control  u(t)  is 
chosen  so  that  it  minimizes  the  quadratic  performance 
index 

J  =  f  f  (xT(t)  Qx(t)  +uT(t)  Ru  (t)  }dt  +  xT(tf)  Sfx(tf)  (2-5) 

where  Q  is  a  time  invariant  positive  semidefinite  state 
penalty  matrix,  R  is  a  time  invariant  positive  definite 
control  penalty  matrix,  and  Sf  =  S(t^) .  This  process 
yields  a  time  invariant  feedback  gain  matrix,  K. 

The  optimal  control,  u*(t),  is  defined  as 

u*(t)  =  -K  x ( t )  (2-6) 

The  feedback  control  gain  matrix  K  is  determined  by  the 
solution  to  the  Riccati  equation: 
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S  =  ATS  +  SA  +  Q  -  SBR~1BTS  S(tf)  =  Sf  (2-7) 

K  =  R~ 1  BT  S  (2-8) 

Therefore,  the  optimal  control  law  is  a  function  of  the 
matrices  A  and  B  defined  in  the  linear  system  model,  and 
of  the  matrices  Q,  R,  and  S ^ ,  defined  in  the  quadratic 
performance  index.  If  the  system  is  both  detectable  and 
stabilizable,  then  as  the  system  reaches  steady  state  (i.e., 
as  t  -*■  00 )  the  solution  to  the  Riccati  Eq  (2-7)  converges 
to  a  constant,  steady  state  value  regardless  of  S^.  This 
solution  is  determined  by  setting  S  =  0.  Thus,  from 
Eq  (2-8)  it  is  clear  that  the  Riccati  feedback  gain  matrix, 
K,  also  becomes  a  constant  and  the  optimal  control  law, 

Eq  (2-6),  becomes  time  invariant  in  steady  state  (Ref  5: 
574-5) .  Therefore,  the  last  term  of  the  quadratic  per¬ 
formance  index,  Eq  (2-5) ,  becomes  negligible  at  large 
enough  values  of  time,  t.  The  term  is  assumed  to  be  zero, 
yielding  the  steady  state  quadratic  performance  index 

Jss  =  fo  {^T(t)  2.  ±{t)  +  u(t)TR  u(t)}  dt  (2-9) 

There  is  no  clearly  defined  method  for  choosing  the 
penalty  matrices  Q  and  R.  Typically,  Q  and  R  are 
initially  chosen  to  be  diagonal  matrices  with  each 
diagonal  element  equal  to  the  reciprocal  of  the  maximum 
allowable  excursion  squared.  Then  Q  and  R  must  be  tuned 
to  achieve  the  desired  response.  The  tuning  procedure 
is  fairly  arbitrary  and  usually  done  by  engineering 
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intuition.  In  general,  the  process  is  an  iterative  one 
with  the  following  steps: 

1.  Choose  values  for  Q  and  R 

2.  Determine  the  Riccati  feedback  gain  matrix,  K 

3.  Determine  the  closed-loop  state  matrix,  (A-BK) 

4.  Determine  the  frequency  response  and/or  time 
response  characteristics  of  the  system 

5.  Repeat  steps  1  through  4  until  the  system  has 
the  desired  response  characteristics 

The  larger  the  value  of  a  particular  term  in  the  Q  matrix, 
the  more  that  particular  state  value  will  be  penalized  for 
deviation  from  some  desired  nominal  trajectory.  Similarly, 
the  larger  the  value  of  the  R  matrix,  the  larger  the 
penalty  for  using  excessive  control. 

Classical  Control  Methods 

Classical  control  methods  typically  use  frequency 
domain  models  and  analysis  techniques  in  the  design  and 
compensation  of  control  systems.  These  methods  use  the 
transfer  function  between  the  input,  u,  and  the  output, 
y,  that  is: 


G (s)  =  y (s)/u(s)  (2-10) 

where  G(s)  is  the  transfer  function  matrix  and 

G . . (s )  =  y.(s)/u.(s).  u,  ( s)  =  0,k/j  (2-11) 

-LJ  1  j  K 

Classical  techniques  typically  include  the  root-locus 
method,  Bode-plot  representations,  Nyquist  diagrams,  and 
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Nichols  charts  to  assist  the  designer  in  giving  a  system 
the  specified  response  characteristics.  Unfortunately , 
these  methods  are  mostly  limited  to  single-input-single¬ 
output  (SISO)  systems,  and  cannot  easily  be  used  for 
multiple-input-multiple-output  (MIMO)  systems.  In  the 
MIMO  system  case,  as  one  transfer  function  is  compensated, 
it  affects  others  and  designing  a  control  system  using 
classical  methods  becomes  increasingly  more  difficult, 
especially  for  coupled  off-diagonal  terms.  Optimal  control 
design  methods,  using  vector  mathematics,  are  much  more 
suited  to  the  MIMO  systems.  Linear  quadratic  regulators 
have  an  advantage  that  the  controller  is  guaranteed  to  be 
stable  with  60°  phase  margin  and  50%  to  infinite  gain 
margin.  A  disadvantage  of  classical  methods  is  that  the 
designer  must  be  careful  to  insure  that  the  system  is 
stable . 

Frequency  Shaping  of  Cost  Functionals 

Gupta  (6:5)  has  suggested  that  optimal  control 
theory  can  also  be  applied  to  cost  functionals  involving 
frequency  dependent  state  and  control  penalty  matrices, 

Q  and  R,  respectively. 

Gupta  states  that  some  of  these  controllers  behave 
like  lag  compensators,  and  that  they  reduce  the  high  fre¬ 
quency  response  of  the  system.  However,  the  frequency 
dependent  penalty  matrices  cannot  be  mixed  outright  with 
the  time  dependent  state  and  control  variables  in  the  per¬ 
formance  index  (Eq  2-9).  Therefore  the  steady  state 
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quadratic  performance  index  is  written  initially  in  the 
frequency  domain  in  the  following  manner: 

J  =  f_x  ix*  ( j  uj)  Q  (oj)  x  ( jw)  +  u*  (jw)R(ui)u(jw)  }  du  (2-12) 

where  the  asterisk,  *,  means  complex  conjugate  transpose, 
and  the  Q(x)  and  R(co)  state  and  control  penalty  matrices, 
respectively,  are  specified  to  be  Hermitian  positive  semi- 
definite  and  Hermitian  positive  definite,  respectively. 

Once  the  frequency  domain  performance  index  is  defined, 
it  is  transformed  to  the  more  readily  usable  time  domain 
performance  index  using  Parseval's  Theorem.  That  is,  the 
performance  index  is  evaluated  by  using  the  time  domain 
performance  index  transformed  from  the  specified  frequency 
domain  performance  index.  The  Riccati  matrix  feedback 
gain  is  calculated,  which,  in  turn,  gives  the  optimal  con¬ 
trol,  u* ,  for  the  closed  loop  system.  The  closed-loop 
plant  matrix,  A^L  =  (A  -  B  K) ,  is  used  to  analyze  the 
closed-loop  performance  of  the  system. 

Control  Law  Design 

Having  defined  the  frequency  domain  performance 
index  (Eq  2-12),  the  next  step  is  to  factor  the  terms 

1* 

x*  ( jw)Q(w)x(  ju)  and  u*(jx)  R(w)u(juj)  into  the  terms  x  (jw) 
1  1*1 

x  (j,0  and  u  (jw)u  ( j respectively  where 

^(jw)  =  P1  (  jco)  x  (  joj)  (2-13) 

u1(j'.)  =  P2(jw)u(jw)  (2-14) 
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I  and 


P1*(juj)  P  x  ( j  uj  )  =  Q(  jw) 


(2-15) 


P2*(jw)  P2(jca)  =  R(jw)  (2-16) 

The  pair  of  matrices  P^*(ju)  and  P  (jw)  and  the  pair 
P2*(j°J)  an<^  P_2  ( j 00 )  are  complex  conjugate  transpose 
matrices.  The  following  relationship  exists  between  x(jw) 
and  x^ ( jw) : 


x  ( juj) 


( jcJ) 


xL(j; 


That  is,  we  get  x^(jw)  by  putting  x(jw)  through  a  linear 
system,  G  (jw) ,  whose  frequency  response  is  described  by 
the  transfer  function  matrix  P^(jw),  i.e., 


— x  ( j  °° )  =  i-i ( ^ w ) 


(2-17) 


This  spectral  factorization  of  Q(w)  and  R(w)  is 
necessary  to  facilitate  the  Fourier  transformation  of 
the  frequency  domain  performance  index  into  the  time 
domain.  That  is, 


Q(w)  =  P1*(ju)  P-^jw)  (2-18) 

R(w)  =  P2*(jw)  P2(ju)  (2-19) 

and 

J  =  { x*  ( j  - )  Pj_*(jw)  P-^Cjw)  x(jw) 

+  u*(jw)  P2*(jw)  P  2  (  j  w )  u(jw)  }  d...  (2-20) 
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Therefore  the  performance  index  becomes 

J  =  ,'_x  {x"*-  (jcu)  I  x 1  ( j  ^ )  +  u1  (jou)  I  u^  ( jx)  }du  (2-21) 
1*  1 

where  x  ,  x  ,  and  the  first  _!  are  n  x  n  matrices  and 
1*1 

u  ,  u  ,  and  the  second  I_  are  m  x  m  matrices.  Using  the 
inverse  Fourier  transformation,  Eq  (2-21)  can  be  written 
equivalently  in  the  time  domain  as  the  following: 

oo  .  1 T  ]  ]_T  T 

J  =  .Q  lx  (t)  I  xx(t)  +  ux  (t)  I  ux ( t) }dt  (2-22) 

Using  the  relation  of  Eq  (2-13),  an  equivalent  model  of 
the  system  can  be  described  by  the  following  differential 
equations  with  output  x^ ,  state  z^ ,  and  input  x. 


z.ji  (t)  =  (t) 

+  B-^xft) 

(2-23) 

xX(t)  =  C^tt) 

+  E^xU) 

(2-24) 

The  feed- 

forward  matrix,  D^,  is 

zero  if  P^ ( j  w) 

is  proper, 

that  is , 

if  the  degree  of  the  P^ 

(jo)  numerator 

is  of 

smaller  order  than  the  P^(ju)  denominator.  Eqs  (2-23)  and 
(2-24)  model  the  augmented  system  which  includes  the  ori¬ 
ginal  states  and  the  frequency  weightings  of  the  con¬ 
troller.  Appendix  D  illustrates  the  use  of  the  matrix. 
A  similar  relation  exists  for  the  control  relation  des¬ 
cribed  in  Eq  (2-14)  with  output  u"'" ,  state  z^,  and  input  u: 


z2  (t) 

=  A2z2(t) 

+  B2u(t) 

(2-25) 

u1  (t) 

=  C2z2(t) 

+  D2u(t) 

(2-26) 
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Fig.  2-2  shows  a  block  diagram  of  the  generalized  con¬ 
troller  (Ref  5:531).  Therefore  the  entire  system  dynamics 
can  be  described  by  the  following: 


d 

dt 


’x(t)“ 

1 

o/ 

ol 

<1 

1 _ 

’x(t) 

-  - 

B 

z_i  (t) 

= 

Sx  0 

!]_ 

+ 

0 

Z-  (t) 

_  ^  - 

_»  »  a  2_ 

2  * 1 1 

2. 

(2-27) 


Eq  (2-27)  can  be  defined  in  terms  of  the  following  new 
variables  and  matrices: 


x(t)  =  [A]x(t)  +  [  B  ]  u  ( t ) 


(2-28) 


The  new  steady  state  performance  index  becomes: 


J„  =  /*{  xT  ( t)  Qx  ( t)  +  uT (t) Ru(t) }dt 
ss  u 


(2-29) 


where  Q  and  R  are  the  state  and  control  penalty  matrices, 
respectively,  formed  from  the  augmented  system  defined  in 
Eqs  (2-23)  through  (2-26).  Eq  (2-29)  can  be  written  in 
matrix  notation  to  be  the  following: 


.co  T  ry 

Jss  =  -'g  [x  (t)  u  (t)  ] 


“  * 

Q  0 

x(t) 

0  R 

u  ( t) 

(2-30) 


or  equivalently: 


T  T  1 
d:d,  D,C,, 


T  ,  T,,  .  T,  ,  T , ,  ,  T,..,  -1-1  -1-1, 

Jss  =  J0  -  ^  '*2^  '-  - 1 


0 

0 


0  0 

■  x  ( t ) 

0  o 

z1  (t) 

T  T 

C  C  CD 

—2—2  -2—2 

T 

i2(t) 

T  T 

DC  D  D 

u(t) 

2  2  2  2-- 

—  — ' 

dt 


(2-31) 
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Solving  Eqs  (2-27)  and  (2-31)  in  the  steady  state  Riccati 
Equation  yelds  S  which  defines  the  Riccati  feedback  gain 
matrix  K  for  the  augmented  system.  The  optimal  control 
for  the  augmented  system  is  defined  to  be: 

u*(t)  =  K  x { t )  =  K  x  +  K  z.  +  K  z_  (2-32) 

"~x  ?  2  x  z  ^ 

The  scope  of  this  report  is  limited  to  the  analysis  of  the 
time  response  and  Bode  frequency  response  characteristics 
of  the  closed-loop  system. 
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Ill  Comouter  Programs  and  Methods 


Once  the  theory  of  incorporating  frequency  weighted 
cost  functionals  is  established,  it  is  necessary  to  imple¬ 
ment  the  closed-loop  augmented  system  in  some  readily 
usable  manner.  The  synthesis  of  the  new  system  is  done 
with  the  use  of  digital  computer  programs.  The  use  of 
these  programs  is  described  in  this  section. 

Use  of  OPTCON  Program 

Having  defined  the  augmented  system  matrices  A,  B, 

Q,  and  R  the  first  step  in  analyzing  the  augmented  system 

is  to  calculate  the  Riccati  matrix  feedback  gain,  K  = 

[K  ,  K  ,  K  ] .  There  are  currently  several  programs 
x  z1  z2 

which  solve  the  steady  state  Riccati  equation.  The  pro¬ 
gram  MRIC  by  Kleinman  (Ref  9)  can  be  used  with  other  sub¬ 
routines  to  multiply,  add,  invert,  and  transpose  matrices 
to  solve  the  Riccati  equation.  The  program  OPTCON  (Ref  14) 
is  an  interactive  program  for  optimal  control  analysis 
that  incorporates  MRIC  and  these  other  matrix  operator 
subroutines  to  solve  the  Riccati  equation.  The  new 
system  matrices  A  and  B  and  the  penalty  matrices  Q  and  R 
are  entered  into  OPTCON.  OPTCON  will  output  the  Riccati 
solution,  Riccati  feedback  gain  matrix,  K,  and  the  closed- 
loop  eigenvalues. 

Use  of  TOTAL  Program 

Once  the  Riccati  feedback  gain  matrix,  K,  is  deter¬ 
mined,  the  new  closed-loop  system  plant  matrix, 


is 


r 


formed  as  A^L  =  (A  -  BK) .  Now  the  Bode  frequency  response 
and  time  response  of  the  closed-loop  system  is  analyzed 
using  the  interactive  program  TOTAL  (Ref  9).  The  A^ 

C,  and  D  matrices  are  input  and  the  closed-loop  transfer 
function  of  the  new  system  is  calculated.  Then  the  Bode 
frequency  response  characteristics  and  the  time  response 
characteristics  of  each  system  can  be  analyzed.  This 
process  is  repeated  for  each  of  the  various  cases  to 
determine  any  advantages  of  one  controller  over  the 
others . 
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IV  Investigation  of  Frequency 
Weighted  Cost  Functionals 

A  simple  open-loop  unaugmented  baseline  system  is 
analyzed.  It  is  then  augmented  with  a  baseline  controller 
whose  penalty  matrices  are  identity  matrices.  The  per¬ 
formance  and  characteristics  of  the  closed-loop  baseline 
system  with  the  baseline  controller  is  compared  to  several 
of  the  frequency  weighted  controllers  proposed  by  Gupta. 
The  two  types  of  controllers  that  are  examined  use  fre¬ 
quency  dependent  state  penalty  matrices  and  frequency 
dependent  control  penalty  matrices.  The  two  controllers 
are  examined  independently  to  isolate  the  effect  of  the 
state  or  control  penalty  matrix  on  the  performance  index. 
The  Bode  frequency  response  and  time  response  characteris¬ 
tics  of  each  closed-loop  system  is  analyzed  and  compared 
with  the  baseline  controller  closed-loop  system. 

Baseline  System 

The  baseline  model  is  a  second  order  linear  system 
with  an  open-loop  natural  frequency  of  =  1.0  and  damp¬ 
ing  of  ;  =  0.2.  This  simple  lightly  damped  system  is 
easy  to  analyze  in  terms  of  natural  frequency  and  changes 
in  damping.  The  corresponding  characteristic  equation 
for  this  system  is  the  following: 

s2  +  0.4s  +  1.0  =  0  (4-1) 

The  state  equations  for  this  system,  expressed  in  phase 
variable  control  canonical  form,  are  the  following: 
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d 

dt 


-x(t)- 

-X  (t). 

X  (  t) 

0 

+ 

X  (  t ) 

1 

u  ( t) 


:  4-2) 


Z(t)  =  [1  0]  x(t) 


(4-3) 


The  transfer  function  between  the  state  x  and  the  control 
u  is  calculated.  Fig.  4-1  shows  the  Bode  plot  for  the 
baseline  system.  The  open-loop  system  has  a  high  fre¬ 
quency  slope  of  -40  dB/decade.  The  baseline  system  time 
response  is  shown  in  Fig.  4-2.  The  system  has  a  peak 

time  of  t  =3.2  sec,  peak  overshoot  of  53%,  settling  time 
P 

t  =  19.6  sec  and  final  value  of  1.00. 
s 

Baseline  System  with  Baseline  Controller 

The  baseline  controller  uses  state  and  control 
penalty  matrices  that  are  identity  matrices.  That  is: 


and 


J 


ss 


Q  = 


1 

0 


O' 

1. 


R  =  [1] 


(4-4) 


-'ll  ( x ' *  (  j  u ) 


x  ( jw) 


+ 


u  (jw)  fl]  u  (jus)  }du)  (4-5) 


While  this  simple  case  with  identity  Q  and  R  matrices  is 
not  necessarily  the  best  conventional  controller,  it 
serves  as  a  good  model  to  be  compared  with  the  controllers 
that  use  frequency  weighted  cost  functionals.  Table  4-1 
shows  the  frequency  response  characteristics  for  the 
baseline  controller. 
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FREQUENCY  RESPONSE  OF  OPEN  LOOP  BASELINE  SYSTEM 


o 


Frequency  Response  of  Open  Loop  Baseline  System 


ME  RESPONSE  OF  OPEN  LOOP  BASELINE  SYSTEM 


Time  Response  of  Open  Loop  Baseline  System 


V 


Table  4- 

1.  Characteristics  of 

Baseline 

Controller 

I 

DC  Gain 

|  Closed-Loop  Transfer 
Function 

Gain  |  Poles 

Natural 
Freq,  wn 

i  ; 

Damping 

Ratio 

r 

V 

0.708 

1 

1.00  -0 . 705^0 . 958 j 

1.19 

0.593 

This  baseline  controller  system  increased  the  damp¬ 
ing  to  ;  =  0.593,  and  had  a  natural  frequency  of  ^  =  1.19. 

The  Bode  plot  for  this  system  is  shown  in  Fig.  4-3.  The 
dc  gain  is  0.708  and  the  high  frequency  rolloff  slope  is 
-40  dB/decade.  Fig.  4-4  shows  the  time  response  of  the 
baseline  controller  system.  Using  this  system,  the  time 
response  performance  of  the  system  was  improved  to  the 
values  shown  in  Table  4-2. 

Table  4-2.  Baseline  Controller  Time  Response  Characteris¬ 
tics 


i 

l  Rise  Time 

Peak  Time 

Settling  Time 

Peak 

I  Final  ! 

i  (sec) 

(sec) 

( sec) 

Value 

,  Value 

1.54 

3.28 

4.99 

0.777 

0.707 

Now  the  use  of  frequency  weighted  cost  functionals  is 
examined . 

Case  A 

The  first  frequency  weighted  cost  functional  models 
a  state  penalty  matrix  that  affects  only  the  state, 
that  is: 
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FREQUENCY  RESPONSE  OF  BRSELINE  CONTROLLER 


o 
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FREQUENCY  (RflO/SEC) 


ME  RESPONSE  OF  BRSELINE  CONTROLLER 


r 


Q(w)  = 


2  2 
w  +  a 


:4-6) 


and 


ss 


—  CO  — 


(x* ( j  u ) 


2  2 
+a 


x(j-)  +  u(j^)  [1]  u(j^)}d^  (4-7) 


Because  this  frequency  weighting  function  has  a  numerator 
of  smaller  degree  than  the  denominator,  it  is  easy  to  form 
the  realization  for  this  system  and  implement  the  controller . 
Also,  this  Q  matrix  only  affects  the  x^  state  and  will  be 
compared  later  to  the  effect  of  frequency  weighting  on 
both  states.  The  Q  matrix  must  be  decomposed  into  the 
product  P^*(jw)  P-^jw)  as  described  by  Eq  (2-15).  This 
is  done  by  the  technique  of  spectral  factorization  (Ref  3: 
507-8) .  The  complex  conjugate  matrix  factors  of  Q  are 
the  following: 


P  x  (  j  U) )  = 


3^+a 

0 


P*  ( jou) 


-  ja)  +  a 
0 


(4-8) 


Therefore,  the  new  variable  x^(jw)  has  been  defined  as 
follows : 


1  ,  .  , 
x  (310)  = 


Recalling  the  augmented  system  equations 

zL(t)  =  A]_z1(t)  +  B1x(t)  (4-10) 


3tj+a 

0 


x  (  j  X  ) 


(4-9) 
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x1  (t)  =  Cj^U)  +  Dxx(t)  (4-11) 

and  if  x^(s)  =  H(s)x(s)  according  to  Eqs  (4-10)  and  (4-11) 
then: 


H  (s )  =  C-^sI  -  Ax)  _1  B1  +  Dx  (4-12) 

However,  instead  of  knowing  the  system  matrices  A^,  ,  C ^ , 

and  and  then  solving  for  the  transfer  function  matrix, 
H(s) ,  the  reverse  is  true.  That  is,  the  transfer  function 
matrix  between  x^(s)  and  x(s)  has  been  defined  by  Eq  (4-9) 
and  the  system  matrices  must  be  solved  for.  A  realization 
(Ref  5:441-4)  for  this  transfer  function  must  be  deter¬ 
mined.  Appendix  B  describes  how  to  create  a  realization 
for  a  system.  The  following  describes  this  particular 
realization: 


ll  <s> 


1 

s+a 

0 


0 

1 


zx  (t)  =  [z_x  (t)  ] 


x  (t) 


xx  (t) 
x2  (t) 


(4-13) 


(4-14) 


£i  ( t)  =  [-a]  z_x  ( t )  +  [1  0]  x  ( t ) 


x1  (t) 


zx (t)  + 


X  ( t ) 


where  A^  =  [-a],  B^  =  (1  0],  =  [1  0]  ,  and  = 

are  the  matrices  described  in  the  new  system  model, 
Eqs  (4-10)  and  (4-11).  Then  we  get  the  following: 


(4-15) 

(4-16) 
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(4-17) 


H { s)  =  C1  (si  -  Ax)  1  B1  +  D1 


a<s>  *  [J  sfe]  11 01 +  B  ?] 


H  ( S )  = 


0 


s+a 
0  1 


(4-18) 


(4-19) 


Eq  (4-19)  checks  with  Eq  (4-9)  to  show  that 


x1  (s) 


(s) 


s+a 
X2  (s) 


(4-20) 


This  new  system  model  can  be  expressed  in  terms  of  both 
the  original  and  new  systems  by  putting  the  model  in  the 
form  of  Eq  (2-27) .  Define  the  augmented  state  vector 
X  ( t )  =  [  x  ( t )  (t) ]T. 


x(t)  =  A  x  ( t)  +  B  u  (t) 


(4-21) 


x  (t)  = 


'a 

0  ' 

'b' 

A 

*1. 

x(t)  + 

0 

U  (  t) 


(4-22) 


y_( t)  =  fc1  o  J  x ( t) 

Therefore,  the  new  system  model  is  the  following: 


x  (t)  = 


1  o  i  !  ol 

r°l 

= 

1-1  -0.4  !  0 

x  ( t )  + 

1 

L  1  0  !  -aJ 

LoJ 

u  (t) 


(4-23) 


(4-24) 


£(t)  =[100]  x(t) 


(4-25) 
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The  state  penalty  matrix,  Q,  is  determined  from  the  steady 
state  performance  Eq  (2-31).  Since  only  the  stare  (and 
not  the  control)  is  being  frequency  shaped  in  this  case, 
only  the  upper  left  quarter  of  the  penalty  matrix  in  Eq 
(2-31)  is  needed.  Thus  the  penalty  matrices  are  the 
following : 


T 

—1—1 

T 

£l£l 

'0 

0 

o' 

Q  = 

= 

0 

1 

0 

.0 

0 

1. 

R  =  [R]  =  [1]  (4-27) 


The  controllability  matrix,  M  ,  for  this  realization  is 
the  following: 


M 


[B,  AB, 


A2B] 


0  1  -0.4 

1  -0.4  -0.84 

.0  0  1  . 


(4-28) 


Since  is  of  full  rank,  the  augmented  system  is  completely 
controllable . 

Eq  (4-24)  represents  the  minimal  realization  of  the 
system.  Fortmann  and  Hitz  describe  a  method  to  form  the 
realization  for  a  given  transfer  function  (Ref  5:441-446). 
Their  method  is  relatively  easy,  but  it  generally  does  not 
yield  the  minimal  realization  for  a  system.  However,  by 
knowing  the  non-minimal  realization,  it  is  easy  to  find 
the  minimal  realization.  The  non-minimal  matrices  are  par¬ 
titioned  to  rank  which  yields  the  minimal  realization. 

The  non-minimal  realization  of  the  system  is  the  following: 
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4 


— 

“ 

0 

1  1 

0 

0 

'O' 

-1 

-0.4 

_0  _ 

0 

x  ( t)  + 

1 

1 

0 

-a 

0 

'o' 

0 

1  i 

0 

-a_ 

.0. 

(4-29) 


±(t)  =  [1000]  x(t) 


(4-30) 


The  difference  between  the  above  realization  and  the  minimal 
realization  is  that  the  Riccati  feedback  gain  matrix,  K, 
is  not  of  full  rank  and  the  system  is  controllable  but 
not  completely  observable. 

The  value  of  "a"  ranged  from  0.1  to  10.0  or  1/10  to 
10  times  the  natural  frequency  of  the  open-loop  baseline 
system.  Table  4-3  shows  the  Riccati  feedback  gain,  closed- 
loop  transfer  function  between  the  state  x,  and  input  u, 
natural  frequency,  ,  damping,  and  dc  gain  for  the 
various  values  of  "a". 

Bode  plots  for  Case  A  are  shown  in  Fig.  4-5.  Each 
of  the  Case  A  augmented  systems  has  an  increased  damping 
ratio.  Compared  with  the  baseline  controller  system,  which 
has  a  damping  ratio  of  ;  r  0.593,  the  augmented  systems  with 
frequency  shapings  have  dampings  that  range  from  0.539  to 
0.589.  Like  the  baseline  system,  each  of  the  augmented 
systems  has  a  high  frequency  asymptote  of  -40  dB/decade. 
However,  with  a  frequency  weighting  value,  "a",  at  a  larger 
value  than  the  natural  frequency,  the  magnitude  response 
percent-overshoot,  in  the  time  domain,  to  a  step  input  is 
much  less.  For  a  =  10  the  frequency  response  is  nearly 
flat  until  the  break  frequency  where  it  drops  off  sharply 
at  -40  dB/decade. 
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FREQUENCY  (RRO/SEC) 


Fig.  4-6  shows  the  time  responses  for  Case  A.  The 
time  response  for  each  of  the  augmented  systems,  with  a 
unit  step  input,  is  much  better  than  the  original  system. 
The  systems  with  tne  larger  values  of  "a"  have  increased 
damping  over  those  with  smaller  values  of  "a"  as  shown  in 
Table  4-4.  The  dc  gains  for  the  different  cases  vary 
greatly  and  result  in  the  large  difference  in  time  response 
final  values. 

Table  4-4.  Time  Response  Characteristics  for  System  A 


>  1 

1  «l  _  11 

d  j 

Rise 
Time 
(sec ) 

Peak 

Time 

(sec) 

Settling 

Time 

(sec) 

Peak 

Value 

Final 

Value 

10 

1.71 

3.72 

5.79 

1.13 

0.995 

5 

1 

1.70 

3.70 

5.75 

1.11 

0.981 

0.5 

0.989 

2.65 

7.25 

0.630 

0.447 

H 

O 

0.347 

2.11 

10.2 

0.442 

0.0995 

1  |  Baseline 

I  (Controller! 

1.54 

3.28 

4 .99 

0.777 

0. 707 

! 

Table  4-4  shows  that  the  closed-loop  augmented  system 
zeros  are  slightly  smaller  than  the  poles.  In  this  regard, 
the  controller  acts  like  a  lead  filter,  but  over  a  very 
small  frequency  band.  Therefore  the  effect  of  the  zero- 
pole  relationship  is  minimized.  In  general,  the  system 
performs  best  both  in  terms  of  frequency  and  time  domain 
characteristics  with  values  of  "a"  that  are  larger  than 
the  open-loop  natural  frequency.  These  results  are  not 
what  I  expected.  For  this  frequency  shaping  of  l/(.  +a  ), 


■  I  ME  RESPONSE  FOR  CASE 


Time  Responses  for  Case 


I  expected  that  as  j  -  0  the  low  frequency  weighting  on 

2  2 
x ^  would  approach  1/a  .  Similarly  as  w  -*■  *  the  ..  term 

would  be  the  overwhelming  term  in  the  denominator  and  that 

the  high  frequency  weighting  would  approach  zero  and  the 

closed-loop  high  frequency  response  would  be  similar  to  the 

open-loop  high  frequency  response.  As  Figs.  4-5  and  4-6 

illustrate,  all  of  the  closed-loop  responses  have  a  larger 

damping  ratio  at  all  frequencies  than  the  open-loop  system. 


Case  B 

The  second  frequency  weighted  cost  functional  models 
a  penalty  function  that  affects  only  the  x^  state,  as  before, 
but  a  numerator  constant  equal  to  the  denominator  constant, 
"a"  is  added  as  described  below: 


Q(w) 


■x  +a 
0 


0 

1 


(4-31) 


This  frequency  shaping  is  identical  to  case  A  with  the 

2 

addition  of  the  a  term  in  the  numerator.  I  hoped  that 
adding  the  numerator  constant  would  cause  the  dc  gain  of 
the  closed-loop  systems  to  be  the  same  for  the  different 
values  of  "a".  Experience  with  classical  frequency  domain 
compensation  techniques  has  shown  that  adding  a  numerator 
constant  equal  to  the  denominator  constant  normalizes  the 
dc  gain  as  desired.  In  the  time  domain  this  has  the  sig¬ 
nificance  of  normalizing  the  final  response  to  a  given 
input.  I  hoped  that  this  Q(x)  matrix  would  make  it 
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easier  to  compare  the  performance  of  the  closed-loop 

system  with  the  different  values  of  "a".  Decomposing 

* 

Q(u>)  into  P^(jw)  and  P^(jw)  we  get  the  following: 


P^jw)  = 


juj+a 

0 


?*  (  joi) 


0 


-j  u,+a 
0  1 


Recalling  the  augmented  system  equations: 


(4-32) 


£1(t)  =  A^Z^i  ( t)  +  B-^^t) 


(4-33) 


x1  ( t)  =  C-^z^t)  +  D^x  (t) 


(4-34) 


a  realization  for  the  system  is  defined  to  be  the 
following : 


z^ft)  =  [-a]  z^U)  +  (a  0]  x(t) 


(4-35) 


xx(t)  = 


J 


z  1  ( t )  .  ^ 


0  0 
0  1 


x  ( t ) 


(4-36) 


Checking  the  transfer  function  between  x  (s)  and  x(s)  des¬ 
cribed  by  Eq  (4-32)  we  see  that  this  is  in  fact  a  realiza¬ 
tion  of  the  system. 


P^  (s)  =  (sl_  ~  h.i 


-  f1  — " 

[0  s+a. 


s+a 
0  0 


1  B1  +  D1 

(4-37) 
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(4-39) 
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T 

Defining  the  augmented  system  state  variable  x  =  (x  z^) 
the  augmented  system  equations  are  the  following: 


&(t)  =  A  x(t)  +  B  u(t) 


°1 

'b* 

ICO 

!> 

»  _ 
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‘  0 
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1  <  O' 

-0.4  ‘  0 
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'O' 
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""  0"  ,  -a 

yr 

u(t) 


(4-41) 


(4-42) 


(4-43) 


The  augmented  system  penalty  matrix,  Q  is  determined  again 
from  the  steady  state  performance  Eq  (2-31)  .  Again  we 
recall  the  following  definition  for  the  augmented  state 
penalty  matrix: 

Q  = 

The  control  penalty  is  the  following: 

R  =  (R)  =  [1]  (4-45) 

This  realization  is  minimal.  Table  4-5  describes  the  fre¬ 
quency  response  characteristics  for  Case  B.  Fig.  4-7  shows 
Bode  plots  for  Case  B  with  difference  values  of  "a".  The 
frequency  response  characteristics  for  the  different 
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(4-44) 


3  5 


FREQUENCY  RESPONSE  OF  CRSE 


37 


FREQUENCY  ( RRD/SEC  ) 


T 


systems  did  not  differ  significantly.  The  natural  fre¬ 
quency  ranged  from  1.19  to  1.002.  In  general,  the  damping 
ratio  increased  from  ;  =  0.537  for  a  =  0.1  to  ;  =  0.591 
for  a  =  10.  All  of  the  systems  had  a  very  flat  frequency 
response  roll  off.  The  controllers  with  large  values  of 
"a"  have  the  smallest  overshoot.  The  rolloff  asymptote  has 
a  slope  of  -40  dB/decade.  All  of  the  systems  have  the  same 
dc  gain  which  is  what  was  desired. 

The  time  response  characteristics  of  Case  B,  to  a 
unity  step  input,  are  shown  in  Fig.  4-8.  Adding  the  con¬ 
stant,  "a",  to  the  numerator  makes  all  final  values  equal 
to  0.707.  Table  4-6  shows  that  the  time  response 
characteristics  are  very  similar  for  all  values  of  "a". 
Using  values  of  "a”  that  were  less  than  the  open-loop 
natural  frequency  resulted  in  greater  overshoot,  while 
larger  "a"  values  had  minimal  overshoot.  The  frequency 
shapings  for  Case  B  performed  somewhat  similarly  to 

Table  4-6.  Time  Response  Characteristics  for  System  B 


1 

II  _  II 
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Rise 

Time 

(sec) 

Peak 

Time 

(sec) 

Settling 
Time 
(sec ) 

Peak 

Value 

Final  | 

Value  1 

1 

;  io 

1.54 

3.27 

4.98 

0.778 

0.707 

5 

1.53 

3.26 

4.98 

0 . 779 

0.707 

1 

1.42 

3.17 

7.00 

0.824 

0.707 

0.5 

1.33 

3.19 

7.63 

0 . 885 

0.705  ' 

0.1 

1.21 

3.49 

Large 

1.05 

0.707 
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standard  linear  quadratic  regulator  designs.  The  frequency 
shapings  with  larger  values  of  "a"  imposed  a  greater 
penalty  on  state  trajectory  deviation  than  smaller  values 
of  "a".  In  the  neighborhood  of  the  open-loop  natural 
frequency  (e.g.,  from  0.1  to  10.0  rad/sec)  the  similarity 
is  closest.  Using  conventional  linear  quadratic  methods, 
for  w  =  10  the  weighting  for  a  =  10  is 


2 

a 


uj 


2 


100 

100+100 


0.5 


and  the  weighting  for  a  =  0.1  is 


(4-46) 


a2  0.01 

2  2  100+0.01 
u)  +a 


0.0001 


(4-47) 


At  co  =  0.1  the  weighting  for  a  =  10  is 


a2  _  100  ~ 

2 ,  2  0.01+100 

co  +a 


(4-48) 


and  the  weighting  for  a  =  0.1  is 


a2  _  0.01 

22  0 .01+0 .01 

co  +a 


0.50 


(4-49) 


The  smaller  weightings  behave  more  like  the  lightly  damped 
open-loop  system  and  the  larger  weightings  should  result 
in  larger  damping  ratios.  Qualitatively,  the  systems 
using  frequency  weighted  cost  functionals  behaved  similarly 
to  systems  with  conventional  cost  functionals.  Using  larger 
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values  of  the  frequency  weighting  results  in  larger  damp¬ 
ing  ratios. 


Case  C 


The  third  frequency  weighted  cost  functional  to  be 
discussed  is  the  same  as  Case  B  except  that  it  affects 
both  states,  that  is: 


QUO 


2  ,  2 
u  +a 


(4-50) 


This  case  analyzes  the  effect  of  weighting  both  the  and 
*2  states  instead  of  only  x^  in  the  previous  case.  The 
resulting  change  in  the  frequency  and  time  response  from 
the  previous  case  is  studied.  Using  spectral  factoriza¬ 
tion,  Q(w)  is  decomposed  into  the  following  values  of 
P*  (  jw)  and  P ^ ( j  w )  : 


P-j_  ( jw) 


a 

ju+a 

0 


0 

a 

jw+a. 


Pj(jw)  = 


-gu'+a 

0 


-go+aj 


(4-51) 


The  controller  system  state  equations  for  the  realization 
of  this  transfer  function  are  the  following: 


(t) 


(t) 
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(t) 


(4-52) 
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X  ( t) 


(4-53) 
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The  verification  that  this  system  is,  in  fact,  a  realiza¬ 
tion  of  the  transfer  function  between  x^  (S)  and  x(S)  is 
shown  in  Appendix  C.  The  new  system  can  be  described  by 
the  following: 


0 

1  l  0 

0 

0 

-1 

O 

rr 

o 

0 

1 

a 

”  “0“T  ~a~ 

o 

x(  t)  + 

'o' 

0 

a  ;  0 

-a 

0 

(4-54) 


The  new  penalty  matrices  Q  and  R  are  defined  by  Eq  (2-31) 
to  be  the  following: 


Q  = 


D?C. 

'0  0  0  0' 

-1-1  -1—1 

0  0  0  0 

_T_ 

_T_ 

0  0  10 

^ij 

.0  0  0  1. 

R  =  [R]  =  [1] 


(4-55) 


This  realization  is  minimal.  Table  4-7  shows  the  frequency 
response  characteristics  for  Case  C.  Values  of  the  value 
"a",  which  are  larger  than  the  natural  frequency  have  much 
higher  damping  ratios  (;  =  0.591)  than  the  original  system. 
Unlike  the  previous  cases,  however,  smaller  values  of  "a" 
reduce  the  system  damping  ratio  to  nearly  the  same  value 
(:  =  0.210)  as  the  baseline  system.  This  is  because  Q(^) 
approaches  zero  as  "a"  approaches  zero  which  results  in 
the  augmented  system  responding  more  like  the  open-loop 
system.  Fig.  4-9  shows  the  Bode  plots  for  Case  C.  The 
Bode  plots  are  significantly  different  because  of  the  dif¬ 
ferences  in  damping  ratios,  discussed  previously.  Using 
a  value  of  a  =  10  caused  the  system  to  have  a  flat  magni¬ 
tude  response  with  little  overshoot.  Conversely,  the 
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systems  with  values  of  "a"  smaller  than  the  natural  fre¬ 
quency  have  large  magnitude  overshoot  at  the  break  fre¬ 
quency.  Table  4-7  shows  that  the  natural  frequencies  are 
fairly  constant  for  the  different  values  of  "a".  The 
high  frequency  asymptote  for  all  of  the  Case  C  systems 
is  -4G  dB/decade. 

These  results  can  be  expected.  By  imposing  fre¬ 
quency  shaping  on  both  states,  the  maximum  value  that  Q(^) 
can  have  is  unity,  so  shapings  with  large  values  of  "a" 
should  be  similar  to  the  previous  case  which  also  has  a 
maximum  weighting  approaching  unity.  However,  for  small 
values  of  "a",  as  "a"  approaches  zero,  both  state  weightings 
approach  zero.  The  biggest  difference  between  this  case 
and  the  previous  case  should  be,  and  is,  in  systems  involv¬ 
ing  smaller  values  of  "a",  which  results  in  smaller  damping 
ratios . 

The  time  response  characteristics  for  Case  C  are 
shown  in  Fig.  4-10.  Table  4-8  shows  that  for  larger  values 
of  "a",  the  systems  have  less  overshoot  with  good  time-to- 
peak-value  properties.  The  time  response  characteristics 
for  systems  with  small  values  of  "a"  have  smaller  damping 
ratios  and  are  not  much  better  than  the  open-loop  response. 
Frequency  weighting  on  both  states  has  a  much  more  varied 
effect  on  both  frequency  and  time  response  characteristics 
than  the  previous  cases.  Changing  the  values  of  "a"  causes 
very  different  responses,  unlike  Case  B  that  uses  frequency 
weighting  on  only  one  state. 
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Table  4-3. 


Time  Response  Characteristics  for  System  C 


"a" 

Rise 

Time 

(sec) 

Peak 

Time 

(sec) 

Settling 
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(sec) 

Peak 

Value 

Final 

Value 

10 

i  1.54 

3.27 

4.97 

0.778 

0.707 

5 

;  1.51 

3.23 

|  4.93 

0.781 

0.707 

0.5 

'  1.10 

2.90 

}  10.2 

1.03 

0.707 

0.1 

|  0.952 

3.10 

|  23.6 

1.40 

0.707 

Case  D 

The  fourth  frequency  weighted  cost  functional  to  be 
discussed  uses  both  numerator  and  denominator  weighting 
values  on  both  state  terms,  that  is: 


Q(oj) 


2u_  2 

w  +a 


w2+b2 


2  2 
u>  +a 

2  , ,  2  , 
w  +b  J 


(4-56) 


Gupta  suggests  this  improper  weighting  form  because 
it  provides  a  wide  variety  of  ways  to  shape  the  response 
of  the  closed-loop  system.  Q(w)  is  factored  into  the  terms 
P* ( jw)  and  P1 ( jw)  in  the  following  manner: 


Px ( jo) 


juj+a 

n 

[■-jui+a 

0 

-  jw+a 

jcu+b 

0 

ju+a 

Pj(jw)  = 

-  joj+b 

0 

jui+b  J 

-jx+b. 

(4-57) 


The  augmented  system  equations  are  put  into  the  form  of 
Eqs  (2-23)  and  (2-24)  and  the  modified  system  equations 
become  the  following: 
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(4-58) 


(4-59) 


Appendix  D  illustrates  how  to  form  a  realization  for  this 
system.  The  augmented  system  can  be  described  by  the 
following  state  equation: 


'  0 

1 

i  0 

o' 

'O' 

-1 

-0.4 

1  0 

0 

x(t)  + 

1 

1 

”  o' 

~  -b“ 

"0 

"o' 

.  0 

1 

l  0 

-b. 

.0. 

The  augmented  state  penalty  matrix  is  defined  by  the  and 
matrices  in  Eq  (4-59)  to  be  the  following: 
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This  realization  is  minimal.  Table  4-9  shows  the  charac¬ 
teristics  for  system  D.  Using  both  numerator  and  denomina¬ 
tor  weightings  provides  a  wider  variety  of  ways  to  shape 
the  baseline  system.  Using  a  parameter  value  of  a  =  10, 
larger  than  the  open-loop  natural  frequency,  and  a  value  of 
b  =  0.1,  smaller  than  the  natural  frequency,  the  damping 
ratio  is  increased  significantly  to  '  =  0.657  and  the  new 
natural  frequency  increases  threefold  to  3.19.  Using 
values  of  a  =  1.0  and  b  =  2.0,  the  damping  ratio  is 
increased  moderately  to  :  =  0.4  and  the  natural  frequency, 


48 


u)n  =  1.032  is  almost  the  same  as  the  baseline  system.  For 
the  case  of  a  =  0.1  and  b  =  2,  the  least  damping,  ;  = 

0.172  occurs.  This  damping  is  significantly  less  than 
the  baseline  system. 

Intermediate  values  of  "a"  and  "b"  give  widely  dif¬ 
ferent  damping  ratios  and  natural  frequencies.  Fig.  4-11 
shows  Bode  plots  of  selected  values  of  "a"  and  "b".  The 
system  with  a  =  2,  b  =  1  shows  a  nearly  flat  magnitude 
response  until  the  break  frequency  where  the  slope  drops 
off  at  -40  dB/decade.  All  of  the  high  frequency  asymptotes 
drop  off  at  -40  dB/decade.  This  particular  case  behaves 
similarly  to  a  lead  filter  in  terms  of  the  pole-zero 
location  in  that  the  closed-loop  system  zero  is  less  nega¬ 
tive  than  the  dominant  pole,  and  there  is  a  large  increase 
in  wn  which  reduces  the  settling  time.  The  closed-loop 
transfer  function  has  a  dc  gain  which  is  less  than  unity. 

To  normalize  each  of  the  Case  D  systems,  so  that  each 
system  had  the  same  dc  gain,  the  frequency  weightings 
could  have  been  multiplied  by  a  constant  equal  to  b/a. 
Therefore,  as  seen  in  Fig.  4-11,  the  low  frequency  magni¬ 
tude  could  be  increased  by  multiplying  by  the  additional 
gain . 

Fig.  4-11  and  Table  4-9  show  varying  degrees  of  the 
lead  filter  effect  for  different  values  of  "a"  and  "b"  . 
However,  as  "a"  is  held  constant  and  "b"  is  decreased, 
the  damping  ratio  decreases.  Fig.  4-12  shows  the  time 
response  for  selected  Case  D  systems.  Table  4-10  shows 
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Table  4-10.  Time  Response  Characteristics  for  System  D 
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the  time  response  characteristics  for  Case  D.  The  systems 
with  smaller  values  of  "a"  and  larger  values  of  "b",  (a  = 
0.1,  b  =  10),  have  smaller  damping  ratios,  ;  =  0.212,  and 
greater  values  for  overshoot  and  settling  time.  The  system 
with  a  =  2,  b  =  0.1,  on  the  other  hand,  has  a  much  smaller 
time  response. 

In  general.  Case  D  systems  provide  much  more  flex¬ 
ibility  in  designing  system  response  than  the  previous 
three  cases.  Wide  ranges  of  frequency  and  time  response 
characteristics  can  be  achieved  depending  on  the  choice  of 
"a"  and  "b".  Let  us  assume  that  a  design  objective  is  to 
reject  disturbances.  The  system  can  be  modeled  by  the 
following  block  diagram: 
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f 


Fig.  4-11  shows  Bode  plots  for  the  frequency  response 

between  the  state,  x. ,  and  commanded  control  u  .  If  u  is 

1  c  c 

assumed  to  be  some  disturbance  instead,  then  the  Bode  plots 
can  be  interpreted  as  the  frequency  response  of  the  trans¬ 
fer  function  between  and  some  distrubance.  Comparing 
the  various  Bode  plots  in  terms  of  distrubance  rejection, 
one  characteristic  that  is  desirable  is  small  dc  gains. 

This  indicates  that  reduced  magnitude  response  occurs  at 
lower  frequencies.  If  the  goal  is  to  reject  low  frequency 
disturbances,  then  the  more  negative  the  low  frequency 
response  magnitude  in  decibels,  the  better.  From  Fig.  4-11, 
the  system  with  a  =  2,  b  =  1  is  better  at  rejecting  zero- 
mean  white  Gaussian  noise.  The  system  with  a  =  10,  b  =  0.1 
is  much  better  than  the  baseline  controller  closed-loop 
system  at  rejecting  low  frequency  noise  because  the  low 
frequency  penalty  on  state  is  much  higher.  Small  "a" 
values  and  large  "b"  values  create  a  lightly  damped,  slower 
responding  system,  while  large  "a"  values  and  small  "b" 
values  create  a  highly  damped,  faster  system. 


Case  E 

Unlike  the  previous  cases,  this  case  examines  fre¬ 
quency  weighting  on  the  control  instead  of  the  state. 

One  model  proposed  by  Gupta  (Ref  7:5)  is  the  following: 


R(  ju) 


—  ■“ 

jw  +  a 
a 


(4-62) 


With  the  previous  state  weightings  it  was  not  possible  to 
use  frequency  shapings  with  a  numerator  of  greater  degree 
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than  the  denominator.  However,  by  including  the  control 
in  the  state  vector  it  is  possible  to  realize  the  control 
weighting  given  by  Eq  (4-62)  .  The  purpose  is  to  see  if 
this  type  of  frequency  shaping  on  the  control  can  change 
the  damping  ratio  and  natural  frequency  of  the  closed-loop 
system.  Recalling  Eqs  (2-25)  and  (2-26),  the  augmented 
system  can  be  described  by  the  following: 

z_2(t)  =  A2z_2(t)  +  B2u(t)  (4-63) 

u1 (t)  =  C2z2(t)  +  D2u(t)  (4-64) 

where  the  new  control  variable,  u^(t),  is  a  linear  combina¬ 
tion  of  the  commanded  control,  u(t),  and  some  frequency 
weighted  function  of  u(t).  That  is,  the  new  control 
variable  is  defined  to  be  the  following: 

u1(s)  =  [^]  u(s)  (4-65) 

or  equivalently: 

au1(t)  =  u ( t )  +  au ( t)  (4-66) 

u(t)  =  au(t)  -  au^(t)  (4-67) 

This  case  of  frequency  weighting  is  the  same  as  including 
control  (actuator)  dynamics  into  the  state  equation.  This 
system  is  easily  modeled  and  a  realization  is  easy  to  find. 
The  augmented  system,  with  the  new  control  variable  u^  is 
defined  by  the  following  state  equation: 
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or,  for  this  specific  case: 
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(4-69; 


where  =  [-a]  and  =  [a].  In  order  to  analyze  the 
frequency  response  of  this  case,  the  transfer  function 
between  the  commanded  input  u(S)  and  the  state  x(S)  must 
be  determined.  This  cannot  be  done  outright  from  Eq  (4-69) 
because  of  the  use  of  the  new  control  variable  uA  The 
entire  system  block  diagram  is  shown  in  Fig.  4-13. 


Fig.  4-13.  Block  Diagram  of  Case  E 
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The  loop  must  be  closed  on  u(s)  to  calculate  the  transfer 
function  between  the  state  x(s)  and  the  commanded  input 


U(s).  Appendix  A  shows  that  the  transfer  function  matrix 
is  the  following: 


H  ( s )  = 


(s+a  -  aK^) 

s [ (s+Q . 4 )  (s+a  -  aK^)  -  aK  2)  +  (s+a  -  aK,  -  aK^) 

- (s+a  -  aK^) 

[ (s+0.4)  (s+a  -  aK^)  -  aK^J  +  (s+a  -  aK^  -  aK,  ) 

(4-70) 


The  identity  matrix  will  be  used  for  the  state  and  control 
penalty  matrices,  that  is: 


Q 


'1  0  O' 
0  10 
0  0  1. 


R  =  [1] 


(4-71) 


Table  4--11  shows  the  frequency  response  characteristics  for 
Case  E.  Large  values  of  "a"  (a  =  10)  result  in  larger 
values  of  the  damping  ratio,  £  =  0.208  for  a  =  0.1.  Bode 
plots  for  this  case  are  shown  in  Fig.  4-14.  Each  system 
has  very  similar,  flat  low  frequency  response  characteris¬ 
tics.  All  of  the  systems  have  improved  magnitude  overshoot 
properties  compared  with  the  open-loop  baseline  system. 

All  Case  E  systems  had  larger  overshoot  than  the  closed- 
loop  baseline  controller.  This  is  because  the  baseline 
weightings  are  unity  while  the  frequency  shaped  weightings 
can  only  approach  unit-,  as  "a"  gets  very  large.  The 
a  =  10  system  has  the  least  overshoot  of  all.  All  Case  E 
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Frequency  Response  of  Case 


systems  have  a  high  frequency  asymptote  slope  of  -40  dB/ 
decade . 


The  time  responses  of  Case  E  systems  are  shown  in 
Fig.  4-15.  Time  response  characteristics  are  shown  in 
Table  4-12.  The  system  with  a  =  0.1  has  the  largest  peak 
value  (1.27)  and  longest  settling  time  (22.9  sec).  The 
fastest  system  with  the  least  overshoot  has  a  parameter 
of  a  =  10  with  a  settling  time  of  7.42  sec  and  a  peak  value 
of  0.684. 


Table  4-12. 

Time 

Response 

Characteristics 

for  System 

i 

i  j 

|  "a" 

Rise 

Time 

(sec) 

'  Peak 
Time 
(sec) 

Settling  | 
Time 
(sec) 

Peak 

Value 

Final 

Value 

10 

1.43 

3.19 

7.42 

0.684 

I 

0.577  t 

2 

1.28 

2.97 

7.33 

0.705 

0.577  : 

1 

1 

1.12 

2.79 

9.32 

0.772 

0.577 

0.5 

0.978 

2.79 

12.6 

0.921 

0 . 577 

0.1 

0.843 

2.95 

22.9 

1.27 

0.577 

This  method  of  frequency  shaping  can  be  useful  in 
including  control  dynamics  in  a  system.  In  terms  of  other 
systems,  Case  E  qualitatively  responds  most  like  Case  C 
which  had  the  reciprocal  frequency  weighting  on  both  states. 

Case  F 

The  frequency  weighted  cost  functional  for  this  case 
attempts  to  drive  the  state  penalty  weighting  to  a  very 
large  value  at  some  specified  frequency.  Gupta  suggests 
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Time  Responses  for  Case 


that  this  can  be  done  with  the  following  cost  functionals 
(Ref  7:5): 


QU) 


0 


,  2  2.2 
(a)  -a  ) 


(4-72) 


T 

Q(o)  is  factored  into  the  following  values  of  P^(jcj)  and 
Px ( ju) 


( jw)  =  P^ ( jo) 


(4-73) 


The  controller  system  state  equations  for  a  minimal  realiza¬ 
tion  of  the  transfer  between  x1  (s)  and  x(s)  are  the 
following : 


ii<t)  = 


Xx(t)  = 


0  1 

1 - 

O 

o 

2  n 

2.x  (t)  + 

2  . 

a  0 

a  0 

1  0 

0  0 

z_1  ( t )  + 

0  0 

0  0 

x(t) 


X  ( t ) 


(4-74) 


(4-75) 


The  augmented  system  equations  are  the  following: 


x(t)  = 


0 

- 1_ 

i  ; 
-0.4  ' 

0 

0 

0 

...o 

x(t)  + 

0 

1 

a” 

"-a” 

o 

0* 

0 

a  ; 

0 

-J 

0_ 

u ( t )  (4-76) 


The  augmented  system  penalty  matrices  are  defined  by  Eq 
(4-43)  to  be  the  following: 


'0  0  0  0' 
0  0  0  0 
0  0  10 
.0  0  0  0. 


R  =  [1] 


(4-77) 


Table  4-13  shows  the  characteristics  for  Case  F.  All  of 
the  systems  are  either  poorly  damped  or  have  very  small 
natural  frequencies  or  both.  Each  system  has  two  oscil¬ 
latory  modes  except  the  system  with  a  =  10. 

Fig.  4-16  shows  the  Bode  response  characteristics  of 
the  system.  As  with  all  of  the  previous  cases,  the  high- 
frequency  response  is  the  same  as  the  open-loop  system.  As 
with  Case  D,  the  most  significant  result  is  that  this  fre¬ 
quency  shaping  can  be  used  to  reject  low  frequency  dis¬ 
turbances.  The  frequency  shaping  only  has  an  effect  on 
disturbances  or  modes  that  are  of  lower  frequency  than  the 
natural  frequency  of  the  closed-loop  system.  This  result 

4 

is  achieved  by  removing  the  "a  "  constant  in  the  numerator 
of  the  state  cost  functionals.  The  Bode  plots  shown  in 

4 

Fig.  4-16  describe  systems  with  the  constant  "a  "  present, 
but  when  the  constant  is  removed,  the  dc  gain  for  the 
closed-loop  system  is  very  low  for  small  values  of  "a". 

The  frequency  response  does  not  do  anything  parti¬ 
cularly  out  of  the  ordinary  near  the  natural  frequency  of 
the  open-loop  system.  Gupta's  premise  that  this  type  of 
weighting  is  good  for  supressing  a  disturbance  at  w  =  a 
is  not  evident.  Rather  the  controller  is  good  at  suppress¬ 
ing  disturbances  at  frequencies  less  than  the  natural 
frequency . 
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FREQUENCY  (RflD/SEC) 


I 


Table  4-14  shows  the  time  response  characteristics 
for  Case  F.  The  systems  with  small  values  of  "a"  initially 
respond  in  the  opposite  direction  of  the  final  value.  The 
low  damping  ratios  for  the  systems  are  shown.  The  final 
value  may  be  changed  from  -0.707  to  0.707  by  multiplying 
the  system  closed-loop  transfer  function  between  x(t)  and 
u(t)  by  minus  one.  The  settling  times  for  this  case  are 
much  larger  than  previous  cases.  Fig.  4-17  shows  the 
time  response  for  Case  F. 


Table  4-14 

.  Time 

Response 

Characteristics 

for  Case  F 

i 

Rise 

Peak 

Settling 

Peak 

Final 

"a" 

Time 

Time 

Time 

Value 

Value 

i 

(sec) 

(sec) 

(sec) 

10 

1.26 

3.12 

7.30  ^ 

0.884 

-0.707 

5 

1.25 

3.30 

7 .47 

0.884 

-0.707 

1 

1.05 

4.46 

14.4 

1.002 

-0.707 

0.5 

1.02 

5.84 

18.43 

1.097 

-0.707 

0.1 

12.9 

32.3 

Very’  Large 

0.723 

-0.707 
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V  Conclusions 


Frequency  shaping  of  quadratic  cost  functionals  gives 
the  control  engineer  additional  flexibility  in  designing 
control  systems,  especially  in  suppressing  low  frequency 
modes  and  disturbances.  Using  appropriate  values  for  the 
choice  of  compensator  values  of  "a"  and  "b" ,  the  damping 
ratios  can  be  significantly  increased  over  the  baseline 
open-loop  system,  thereby  decreasing  both  overshoot  and 
settling  time.  The  baseline  conventional  controller  has 
a  larger  value  for  the  damping  ratio  than  the  frequency 
shaped  controllers.  This  is  because  the  latter  has 
weighting  matrix  values  that  only  approached  the  unity 
values  of  the  baseline  controller.  However,  the  frequency 
shaped  controllers  demonstrate  the  ability  to  suppress 
rejections  and  system  modes  which  occur  at  frequencies 
less  than  the  natural  frequency  of  the  open-loop  system. 

The  rise  time  can  also  be  reduced  with  the  frequency 
shaped  controllers.  High  frequency  asymptotes  had  slopes 
of  -40  db/decade  for  the  baseline  system  and  all  of  the 
controllers . 

The  controller  with  both  numerator  and  denominator 
frequency  weighting  {s+a/s+lo)  demonstrated  the  most 
flexibility  in  achieving  a  variety  of  design  character¬ 
istics.  State  penalty  frequency  weighting  can  be  inter¬ 
preted  as  processing  the  state  through  a  linear  shaping 
filter.  Control  penalty  frequency  weighting  is  identical 
to  incorporating  control  actuator  dynamics  into  the  system. 


68 


VI  Recommendations 


Optimal  control  using  frequency  weighted  cost  func¬ 
tionals  provides  the  designer  with  a  wide  spectrum  of 
compensation  characteristics.  This  method  can  improve 
the  response  of  many  control  systems  of  interest  to  the 
Air  Force,  particularly  in  low  frequency  disturbance  sup¬ 
pression.  The  effect  of  frequency  shaping  on  the  sup¬ 
pression  of  low  frequency  disturbances  and  system  modes 
should  be  examined  further. 
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Appendix  A 


Calculation  of  Transfer  Function  Between 

Xi  (s)  and  u  (s)  for  Case  E 
—l  — c 


x(t)  =  Ax ( t)  +  Bu ( t)  +  Bu  (t) 


(A- 1 ) 

(A-2) 


u(s)  =  K-^x  (s)  +  K2u(s) 


u(s)  =  -au(s)  +  aK^x(s)  +  aK^u(s! 


(sI-A)x(s)  =  Bu(s)  +  Bu  (c) 


(s+a  -  aK^)u(s)  =  aK^x(s) 


u(s) 


aK-^x  ( s ) 
(s+a  -  aK. 


(A-3) 
( A-4 ) 
( A-  5 ) 
( A-6 ) 

(A- 7 ) 


B  aK1x(sl 

(si-A)x(s)  -  +  B  ^(s: 


(A-8 ) 


Ba  K- 


(sI-A)  - 


(s+a  -  aK^) 


x  ( s )  =  B  u  (si 


(A- 9) 


x  (s)  = 


H  (s)  = 


(sI-A)  - 


Ba  K1 


(s+a  -  aK- 


-1 


B  u  ( s ) 


x  (s) 
u_  (s) 


( sI-A)  - 


Ba  K 


(s+a  -  aK. 


( A-l 0 ) 


(A- 1 1 ) 


here 


0  1 

o" 

_-l  -0. 

4_'  -  = 

1 

-Cal 


-1  =  tKl  K2]  K3  =  k3 


(A- 1 2 ) 

( A-l 3 ) 


Substituting  in  the  values  of  the  matrices  we  get  the 
following : 
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Appendix  B  Creating  a  Realization  for  a  System 

Fortman  and  Hitz  discuss  the  procedure  of  generating 
a  realization  in  detail  (Ref  4:441-451).  Briefly,  the 
goal  is  to  find  the  system  equations 


X(t) 

=  Ax  ( t) 

+  Bu  ( t) 

(B-l ) 

y  (t) 

=  Cx  ( t ) 

+  Du ( t ) 

(B-2 ) 

that  realize  a  system 

transfer 

function  matrix,  H(s) 

where 

H  ( s ) 

=  C(sl-A) 

-1  B  +  D 

(B-  3 ) 

In  general,  the  realization  will  not  be  unique.  If  the 
state,  x(s) ,  the  control  u(s) ,  and  the  output,  ^(s)  have 
dimensions  n,  m,  and  p  respectively,  then  H(s)  will  be 
a  p  x  m  transfer  function  matrix.  The  first  step  in 
creating  a  realization  is  to  express  the  transfer  function 
matrix  in  a  ratio  of  polynomials  in  s,  that  is: 


H  (s) 


3°  +  Si*  +  02s2  +  ...  Qn.!  B-1 

p  +  p. s  +  p~s2  +  . . .  p  "  sn_l  +  s h 
ro  rl  c2  rn-l 


Note  that  Eq  ( B  —  4 )  is  a  proper  transfer  function.  The 
degree  of  the  numerator  is  less  than  the  denominator. 
For  cases  involving  improper  transfer  functions,  the 
polynomial  equation  must  include  a  feed  forward  matrix, 
D. 
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H  { s)  = 


0  4-  Q  (s)  +  ...  Q  Sn_1 

— o  —1  — n-i 

n_^ 

po  +  pi(s)  +  •**  pn-l  s  +  s 


(B-5) 


Appendix  D  illustrates  the  use  of  the  feed  forward  matrix 
in  improper  transfer  functions. 

The  next  step  is  to  describe  the  system  in  control 
canonical  form  using  the  and  values  from  Eq  (B-4) 
or  (B-5)  . 


x  = 


Om 

Im 

Om  ... 

Om 

Om 

Om 

Im  ... 

Om 

Om 

Om 

di 

Im 

Ml 

0 

-p1Im  - 

p^Im  ... 

-p  .  Im 
'n-l— 

£  =  [0 
*-  — o 

— 1  —2  ' ' 

•  fin-11; 

Om 

Om 

x  + 

• 

Om 

Im 

u  (B-6) 


(B-7 ) 


The  submatrices  are  of  dimension  m  x  m,  denoted  by  the  m 
subscript.  While  this  method  will  provide  a  realization, 
in  general,  it  will  not  be  the  minimal  realization.  In 
forming  the  plant  matrix  in  Eq  (B-6) ,  the  designer  should 
fill  in  the  lower  left  element  first.  Then  the  rest  of 
the  matrix  elements  should  be  filled  as  necessary  to 
achieve  the  correct  matrix  dimension. 

The  final  step  is  to  find  the  minimal  realization. 

A  minimal  system  is  both  completely  controllable  and 
observable.  The  minimal  realization  can  generally  be 
found  by  starting  in  the  upper  left  corner  of  the  A  matrix 
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Appendix  C  Formation  of  Realization  and 
Transfer  Function  for  Case  C 


3<jj+a 

0 


jw+a 


H(s)  = 


Q0  +  QiS  +  ...  Qn_x  s 


P0  +  Pj_s  + 


n-1 


n-1  n 

Pn-l  3 


H(s) 


a  +  (1)  s 


(C-l) 


(C-2) 


(C-3) 


-o  [;g  Po  a 

(C-4) 

z^ft)  =  A^Zjift)  +  E^xft) 

(C-5) 

x^Ct)  =  C^z^(t)  +  D^x(t) 

(C-6 ) 

z1(t)  =  [-pQI2]  z^ft)  +  1 1 2  ]  x(t) 

(C-l) 

x1  ( t:)  =  [Qq  ]  Z±(t) 

( C— 8 ) 

h{t)  ■  [1  -|ix(t)  +  [o  l]s(t) 

(C-9 ) 

i1'1’  ■  [o  a]  il111  +  1  S]s(t) 

(C-10) 

P2  (s)  =  C-^sI  -  A1)~1  B1  +  Dx 

(C-ll) 

^  ■  [o  a]  [T  s^a]-1  [J  D  *  [»  »] 

(C-12) 
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Px  (8) 


1 

(s+a) 2 


^(S)  = 


0 

1 

s+a-1 


Pits) 


0 

1 

s+a 


(C-l 3 ) 


(C-14 ) 


(C-l 5 ) 


There  are  several  possible  realizations  for  a  system 
given  some  transfer  function.  The  realization  described 
by  Eqs  ( C— 9 )  and  (C-10)  is  not  exactly  the  same  as  the 
realization  which  is  used  for  Case  C.  Both  realizations 
define  ^(jui),  given  by  Eq  (C-l). 
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Appendix  D  Formation  of  Realization  and 
Transfer  Function  for  Case  D 


o 

jco+b 

o  i-+— 

ju+b 


(D-l) 


n-1 


PL(s)  = 


— Q  +  Qj(s)  +  fin-!  5 

P0  +  Px(s)  +  •••  +  pn-1  Sn  1  +  sn 


+  E  (D-2) 


P-J^S) 

pa  ol  s+a  oH 

L°  aJ  .  Ti  ol  L 0  s+aJ  ,  fi  d 

-  b  ♦  u)«  L°  u  s  +  b  L°  tl 

( D— 3 ) 

1 _ 1 

II 

o 

ol 

_0  a]  P0  *  b 

(D-4) 

Aitt) 

=  [-pQI2]  ( t)  +  t I2 ]  x(t) 

(D-5) 

x1  (t) 

=  t Qq ^  -1  ^ ^  +  [o  I 

(D-6) 

h  = 

fe-lK-]5 

(D-l) 

1 

X  = 

&b 

( D—  8 ) 

P  !  ( s ) 

=  CitsI  -  Ai)-1  Bx  +  Dx 

(D-9) 

^(s) 

[""a-fa  0  s+b  ol-l  [1  ol  pL  0~| 

[_  0  a-b  0  s+bj  [_0  lj  [_0  1_| 

( D— 1 0 ) 

Pits)  = 


a-b 
0  a- 


0~]  sib  °  [l  ol  fl  o“l 

-bJ  o  -4  L°  Lo  U 

s+b 


( D-l 1 ) 
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